Mutagenesis and Metallic DNA 
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Positions of protons in DNA hydrogen bonds are fundamental for the fidelity of the replication 
process. Because the hydrogen bond is partially covalent the electronic structure of DNA plays an 
important role in the replication process. Current studies of electron transfer in DNA have shown 
that simple pictures based on independent electron physics are not sufficient to account for the 
experimental data. In this paper we propose a model for electron-proton interaction in DNA and 
its relationship to mutagenesis. We argue that mutations will be more frequent in DNA sequences 
where electron derealization occurs. To test our theory we calculate the heat capacity of DNA 
at low temperatures and compare with the available experimental data in hydrated DNA. We find 
that our model is in good agreement with the data. We also propose new measurements of infrared 
absorption which would further test our ideas. 



PACS numbers: 87.16.U, 87.14.G, 87.15.M 

Watson and Crick Q pointed out the importance of 
proton transfer as a possible cause of tautomeric base 
pairs in DNA and its deep relationship to mutagenesis. 
Errors in the replication process of DNA can lead to the 
occurrence of mutations related to genetic diseases and 
some forms of cancer. Lowdin || argued that proton 
transfer is a quantum mechanical effect which has purely 
electrostatic origin. Since then many researchers have 
calculated the potential energy profiles for proton trans- 
fer reactions in isolated base pairs || . These studies have 
shown the importance of the electronic configuration at 
each base pair as a stabilizing factor of the hydrogen 
bond (H-bond). Indeed, recent experiments [Q show the 
occurrence of proton transfer in excited base pairs. The 
problem of electron transfer in DNA has been a source 
of intense debate jf| . Many different experimental stud- 
ies of electron transfer by chemical reaction have shown 
that long range electron transfer in DNA may be rather 
different from the same process in proteins, where single 
electron processes are dominant |j] . In this paper we pro- 
pose a model which combines two essential ingredients in 
the physics of DNA: the electronic degrees of freedom, 
and the protons in the H-bonds along the double helix. 
Our model is different from other proposed models for 
electron-proton interactions in DNA [Q, since we take 
into account («) the quantum mechanical character of 
electron motion, (ii) the charging energy due to electron 
transfer between different base pairs, and (Hi) the speci- 
ficity of the electron-proton processes due to differences 
in base pair sequences. 

Let us consider a periodic DNA sequence. It is well- 
known that the electrons in such DNA can move by the 
overlap of ir orbitals along the DNA helix j^] . In this case 
DNA can be viewed as an one-dimensional (ID) solid. 
Because a periodic DNA provides a periodic potential for 
electron motion, basic quantum mechanics tells us that 
the quantum states of the system are Bloch waves with 
wavenumber k and band index n. To each value of k we 



associate a wavelength A of the electron (k — 2n/\), and 
for each band we associate a molecular or atomic state. 
Band structure calculations |9) indicate that such a DNA 
sequence should be an insulator with a band gap of order 
of a few electron volts. In this case the low energy exci- 
tations of the DNA are excitons. Recent experiments in 
Poly(G)-Poly(C) DNA seem to confirm this picture @. 
However, electronic states in biologically relevant DNA 
(non-periodic) depend considerably on the environmen- 
tal conditions (such as pH, hydration, etc). Although 
the study of periodic DNA sequences is very important, 
it probably does not provide final information on the be- 
havior of electrons in DNA under conditions where mu- 
tagenesis is relevant. 

Recent experiments on long range electron transfer in 
DNA indicate that the introduction of Guaninc-Citosine 
(G-C) base pairs enhance electron transfer O]. This 
situation is analogous to that in semiconductors where 
the introduction of electron donors (acceptors) increases 
the charge conduction by the creation of localized states 
close to the band edges. When a sufficiently large num- 
ber of donors or acceptors is introduced into the DNA 
molecule an impurity band can be formed and conduc- 
tion is possible. This picture seems to be supported by 
recent measurements of DNA resistivity j^]. Another 
possible source of ID conduction is via non-linear exci- 
tations such as solitons j^J. As is well-known, disorder 
in ID systems leads to localization of electron wavefunc- 
tions Ej. However, when the localization length is large 
compared with typical electronic scales, electron transfer 
occurs. Moreover, DNA is aperiodic rather than random, 
and Anderson localization is a weaker effect |lq] . 

We assume from the outset that the charge carriers 
(electrons or solitons) of DNA are delocalized over a re- 
gion of size L such that the energy level spacing Hvf,i/L 
(where vf,i is the Fermi velocity in the strand I) is much 
smaller than the thermal energy ksT (where ks is the 
Boltzmann constant and T is the temperature) . The low 
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energy excitations in this case are particle-hole pairs close 
to the Fermi points at ±kp,i- This simple scenario, how- 
ever, is not enough to describe charge propagation in 
DNA because the charging energies for electron transfer 
between different base pairs is of order of electron volts 
as it is common in organic systems. From this point of 
view DNA is very close to being a Mott insulator instead 
of an ordinary band insulator fll4j. In ID, such a conduc- 
tor is called a Luttinger liquid |Tq]- In a Luttinger liquid 
the low lying excitations are bosonic collective modes of 
spin and charge which have solitonic character. 

Long ago Pauling |0 proposed that H-bonds have im- 
portant covalent character, which has subsequently been 
confirmed [H . Proton states at the H-bonds are very im- 
portant because they stabilize the double helix. Protons 
are much heavier than the electrons and therefore much 
slower. In A-T base pairs there are 2 H-bonds anchoring 
the two strands of DNA, while in the G-C case there are 
3 H-bonds. In principle each proton in the H-bond can 
occupy one of its sides, and therefore there are 4 (8) possi- 
ble states for the protons in A-T (C-G) base pairs. These 
states do not have the same energy because of the elec- 
trostatic interaction between the atoms in the base pairs 
with the protons at the H-bonds. Moreover, the differ- 
ence in energy between the states is usually much larger 
than the thermal energy and therefore we only need to 
consider the ground state and the first excited state. We 
call the excitation energy between these two states e m , 
where m labels the position of the base pair along the 
DNA backbone. In this case the proton Hamiltonian can 
be simply written as 

m 

where a z is a Pauli matrix. In this representation the 
ground state of the proton problem is the pseudo-spin 
state | ft) while the excited state is | JJ-). Notice that e m 
can have only two values corresponding to a A-T or G-C 
base pair. Thus, strictly speaking e m is the energy re- 
quired to produce a tautomeric state on an isolated base 
pair. Because of the large proton mass we disregard pro- 
ton tunneling and treat the proton as a classical particle. 
Notice that at temperature T the probability of having 
a local proton configuration in a mutated state is given 
by the Boltzmann distribution: 

P m (T) cx exp [-e m /(fc B T)] (2) 

corresponding to a thermal jump over the energy barrier. 
At room temperature and with e m of the order of tenths 
of eV this probability is very small and mutation is a rare 
event. 

Let us now consider the coupling between protons and 
7r electrons. This coupling depends strongly on the de- 
grees of freedom of the H-bonds and the ir orbitals. On 
the one hand, the protons in the H-bond have their 



own magnetic spin degrees of freedom and non-magnetic 
pseudo-spin variables associated with their position in 
the H-bond. On the other hand, electron states in DNA 
can be labeled by their momentum hk, their true spin 
magnetic quantum number, and by their two possible 
strand numbers: I = ±1. While the strand and the spin 
quantum numbers of the electrons couple directly to the 
pseudo-spin of the proton, the proton spin couples to the 
electron spin via dipolar forces or contact interactions. 
If we view the single H-bond problem as an impurity in 
a chain (or ladder) then the DNA problem is of the so- 
called multi-channel Kondo type jl9j . Since in this paper 
we are only interested with the stability of the H-bond 
due to the presence of electrons and not with the mag- 
netic properties of DNA, we neglect completely the dipo- 
lar interactions which can be studied separately. More- 
over, the symmetry of the coupling between pseudo-spins 
and electrons is reduced due to the strong anisotropy of 
the proton positions. Therefore, in the simplest model, 
the electron-proton problem reduces to a local coupling 
of value A mj ; between the electronic charge density and 
the H-bond at the base pair m on the strand I: 

in .1 

where p m l is the average electron density of the mth base 
pair in the Ith strand and p m j the actual electronic den- 
sity. The physics of this problem is rather simple. In 
the absence of electronic coupling the H-bond of a given 
base pair m will find itself in the lowest energy state 
which is stable at room temperature. Because the pro- 
ton is coupled to the electrons it can generate electronic 
displacements along the helix with a gain of electrostatic 
energy. This gain of electrostatic energy, in turn, can sta- 
bilize a tautomeric configuration of the protons leading 
to mutation. The main question is whether the gain in 
energy due to the electron-proton interaction is enough 
to compensate for the loss of electrostatic energy between 
the proton and the base pair. Moreover, the coupling be- 
tween H-bonds and electrons also generates an exchange 
interaction between H-bonds, which oscillates in space 
with a wavelength of order of the distance a between 
base pairs. 

The low energy excitations of this Luttinger liquid 
can be studied by the bosonization technique in which 
solitonic modes are bosons with momentum hk, energy 
KvF,i\k\ and strand number I = ±1. In Hamiltonian form 
it can be written as jl6| 

h l = ^ ftvF,i\k\bljb k: i 

k,l 

where bk y i (b\ ,) is the boson annihilation (creation) oper- 
ator and vp.i = VF.iy/l + 2Ui/(ttvf,i) is the renormalized 
Fermi velocity due to the electron-electron interaction 
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strength Ui- In our case the spin degrees of freedom are 
completely decoupled (unless we introduce a direct cou- 
pling between the proton and the electron spin) and do 
not appear explicitly into the problem. In particular, the 
density operator which appears in (p3) can be written as 



Pm,l 



2n\k\ 



Ofc,;e 



where x m is the position of the mth pair along the he- 
lix and e ' = [1 + Ui/(ttvf .z)] 1 ^ 4 - The main gain in us- 
ing the bosonization procedure for this problem is that 
it maps the DNA problem given by the Hamiltonian 
H = Hp + Hl + H e - P into an exactly solvable problem 
in many-body physics called the independent boson model 
po| . Our model describes the relaxation of the Luttinger 
liquid along the helix due to the proton configuration. 

The first result of the model is that the excitation en- 
ergy e m of the H-bond is reduced by an amount A m , 
where 



E 



A 



irv F i + Ui 



(4) 



Thus, in the presence of a local metallic electron fluid one 
gets a reduction of the "natural" bias in the system and 
an increase in the probability (J2J) of finding the system in 
a tautomeric state. It is clear from ([I]) that this can only 
happen when coupling between H-bonds and electrons is 
large or when the Coulomb interaction between electrons 
is weak. In metallic DNA the electron-electron interac- 
tions will be screened by the electron motion and there- 
fore A m can be large enough to induce mutation at room 
temperature. This mechanism for mutation seems to be 
consistent with recent findings relating the existence of 
repetitive sequences in DNA (which tend to delocalize 
the electrons) and mutagenesis associated with evolution 

13- 

In order to test our ideas we have calculated the heat 
capacity of DNA using the above described model. There 
are three main contributors to the DNA specific heat at 
low temperatures: the phonons which give a contribu- 
tion growing like T 3 , the electrons which contribute with 
a term proportional to T, and the protons which con- 
tribute with a Schottky term. The total specific heat 
can be calculated once the DNA sequence is known. 
We assume that the renormalized bias e m = e m — A m 
can be written as a series expansion of the temperature: 
e m w A m + B m T. The temperature dependence of the 
bias is due to the exchange interaction between the H- 
bonds propagated by the electron fluid. In a sample of 
DNA the values of A m and B m are distributed around 
an average value which is given by the average number 
of G-C and A-T pairs. For simplicity we assume here a 
gaussian distribution for both quantities. In Fig. [I] we 
compare the results of our model with the experiments 



reported in Ref . |22|] . The agreement between our model 
and the experiments is good. However, there are serious 
questions about sample quality in experiments dealing 
with DNA and therefore we have also calculated, as a 
further check, the infrared (IR) absorption for DNA us- 
ing many-body techniques. The correlation function of 
interest is the proton spectral function which is defined 
as S m {uj) = — 2Im[G m (u>)] where G m {us) is the Green's 
function of the proton at base pair m as a function of the 
frequency lo. We found that IR absorption at a certain 
site m is controlled by a parameter g m defined as 
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nhvp 



(5) 



At zero temperature and g m < 1, DNA shows an X-ray 
edge effect with a sharp absorption edge. This effect is 
a result of the many-body character of DNA electronic 
conduction. If g m > 1 or the temperature is finite the 
IR absorption becomes smooth. For a specific DNA se- 
quence one has to calculate the absorption of the DNA 
as a whole. In Fig. |^ we show the absorption spectrum 
for one of the DNA sequences of Fig. [l] at finite tem- 
peratures. Observe that no sharp X-ray edge is visible 
because of the superposition of many X-ray edges coming 
from individual protons. There is, however, a suppression 
of the absorption at low frequencies with the formation 
of a "pseudo-gap" . 

In conclusion, we have proposed a model for electron- 
proton coupling and mutagenesis in DNA based on the 
idea of charge derealization. We have shown that the po- 
sition of the protons in DNA can change, leading to muta- 
tions when the ir electrons delocalize. We have compared 
our results for the specific heat at low temperatures with 
the available experimental data and found near quanti- 
tative agreement. We have also proposed new IR absorp- 
tion experiments which can further test our model. 
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ment of Energy. 



[1] Watson, J. D. H. & Crick, F. H. C, Nature 171, 737-738 
(1953). 

[2] Lowdin, P.-O., Adv. Quantum Chem. 2, 213-361 (1965). 
[3] Douhal, A. et al, J. M., Chem. Phys. Lett. 256, 370- 376 
(1996). 

[4] Douhal, A., Kim, S. K. & Zewail, A. H., Nature 378, 260- 
263 (1995); Folmer, D. E. et al, A. W., Chem. Phys. Lett. 
287, 1-7 (1998). 

[5] Wilson, E. K. , C&EN 23, 43 (1999). 

[6] Nunez, M. E., Hall, D. B. & Barton, J. K., Chemistry 
and Biology 6, 85 (1999). 

[7] Steenken, S., Biol. Chem. 378, 1293-1297 (1997). 



3 



[8] Eley, D. D. & Spivey, D. I., Trans. Faraday Soc. 58, 411 
(1962). 

[9] Lewis, J. P. et al, Phys. Rev. B 55, 6880 (1997). 
[10] D. Porath et al, Nature 403, 635 (2000). 
[11] Lewis, F. D. et al, Science 277, 673 (1997), Meggers, E. 

et al, J. Am. Chem. Soc. 120, 12950 (1998). 
[12] Fink, H.-W. & Schonengerg, C, Nature 398, 407 (1999). 
[13] A. H. Castro Neto and A. O. Caldeira, Phys. Rev. E 48, 
4037 (1993). 

[14] Mott, N. F. in Metal-Insulator transition (Taylor & Fran- 
cis, London 1974). 

[15] Vida, J. et at, Phys. Rev. Lett. 83, 3908 (1999). 

[16] Voit, J., Rep. Prog. Phys. 58, 977 (1995). 

[17] Pauling, L. in The nature of the chemical bond (Cornell 
Press, Ithaca 1960). 

[18] Isaacs, E. D. et al, Phys. Rev. Lett. 82, 660 (1999). 

[19] Nozieres, P. & Blandin, A., J. Physique 41, 193 (1980). 

[20] G. Mahan in Many-Particle Physics (Plenum Press, New 
York 1990) pg. 285. 

[21] Moxon, E. R. & Wills, C, Scientific American 280, 94 
(1999). 

[22] Mrevlishvili, G. M. et al, Biophysics 40, 518 (1995). 

FIG. 1. Specific heat of DNA : 1. DNA sample with water 
content no — — 2 M H2O/MBP. The results are presented 
by the continuous curve (theoretical) and o (experimental). 2. 
DNA sample with water content n s = 10 - 12 M H2O/MBP 
(mostly A-form DNA). The results are presented by the dot- 
ted curve (theoretical) and □ (experimental). 3. DNA sam- 
ple with water content n E = 22 - 23 M H2O/MBP (mostly 
B-form DNA). The results are presented by the dashed curve 
(theoretical) and A (experimental). 

FIG. 2. Predicted infrared absorption at finite tempera- 
tures, hvF/a = 2eV . 
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